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Nearly all glass forming liquids display secondary relaxations, dynamical modes seemingly dis- 
tinct from the primary alpha relaxations. We show that accounting for driving force fluctuations 
and the diversity of reconfiguring shapes in the random first order transition theory yields a low 
free energy tail on the activation barrier distribution which shares many of the features ascribed 
to secondary relaxations. While primary relaxation takes place through activated events involv- 
ing compact regions, secondary relaxation corresponding to the tail is governed by more ramified, 
string-like, or percolation-like clusters of particles. These secondary relaxations merge with the 
primary relaxation peak becoming dominant near the dynamical crossover temperature T c , where 
they smooth the transition between continuous dynamics described by mode-coupling theory and 
activated events. 



Diversity, a key feature of glassy systems, is most ap- 
parent in their relaxation properties. Dielectric, mechan- 
ical and calorimetric responses of supercooled liquids are 
not single exponentials in time, but manifest a distri- 
bution of relaxation times. The typical relaxation time 
grows upon cooling the liquid until it exceeds the prepa- 
ration time, yielding a non-equilibrium glass, which can 
still relax but in an age dependent fashion. In addition 
to the main relaxations that are responsible for the glass 
transition, supercooled liquids and structural glasses ex- 
hibit faster motions, some distinct enough in time scale 
from the typical relaxation to be called "secondary" re- 
laxation processes [TJ [2} El SI ES] - These faster motions 
account for only a fraction of the relaxation amplitude 
in the liquid but become dominant features in the relax- 
ation of otherwise frozen glass, where they are important 
to the mechanical properties. These secondary relaxation 
processes in the solvation shell of proteins are also promi- 
nent in protein dynamics [BJ. 

The phenomenology of secondary relaxation has been 
much discussed but, owing especially to the problem of 
how to subtract the main peak, the patterns observed 
seem to be more complex and system specific than those 
for the main glassy relaxation. Some of the secondary 
relaxation motions are, doubtless, chemically specific, 
occurring on the shortest length scales. Nevertheless 
the presence of secondary relaxation in glassy systems 
is nearly universal [7J. In this paper we will show how 
secondary relaxations naturally arise in the random first 
order transition (RFOT) theory of glasses and are pre- 
dicted to scale in intensity and frequency in a manner 
consistent with observation. 

The RFOT theory is based on the notion that there 
is a diversity of locally frozen free energy minima that 



can inter-convert via activated transitions. The inter- 
conversions are driven by an extensive configurational 
entropy. RFOT theory accounts for the well known cor- 
relations between the primary relaxation time scale in 
supercooled liquids and thermodynamics 110] as well 
as the aging behavior in the glassy statepT]. By taking 
account of local fluctuations in the driving force, RFOT 
theory also gives a good account of the breadth of the 
rate distribution of the main relaxation US]- Here 
we will argue that RFOT theory suggests, universally, 
a secondary relaxation also will appear and that its in- 
tensity and shape depends on the configurational ther- 
modynamics of the liquid. This relaxation corresponds 
with the low free energy tail of the activation barrier 
distribution. The distinct character of this tail comes 
about because the geometry of the reconfiguring regions 
for low barrier transitions is different from that of those 
rearranging regions responsible for the main relaxation. 
Near to the laboratory T 9 , the primary relaxation pro- 
cess involves reconfiguring a rather compact cluster, but 
the reconfiguring clusters become more ramified as the 
temperature is raised and eventually resembling perco- 
lation clusters or strings near the dynamical crossover 
to mode coupling behavior, identified with the onset of 
non-activated motions[T3]. Reconfiguration events of the 
more extended type are more susceptible to fluctuations 
in the local driving force, even away from the crossover. 
These ramified or "stringy" reconfiguration events thus 
dominate the low barrier tail of the activation energy 
distribution. 

When the shape distribution of reconfiguration pro- 
cesses is accounted for, a simple statistical computation 
shows that a two peaked distribution of barriers can arise. 
This calculation motivates a more explicit but approxi- 
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mate theory that gives analytical expressions for the dis- 
tribution of relaxation times in the tail. In keeping with 
experiment, the theory predicts the secondary relaxation 
motions are actually most numerous near the crossover, 
but of course, merge in frequency with the main relax- 
ation peak in time scale also at that crossover. Further- 
more the relaxation time distribution for secondary re- 
laxations is predicted to be described by an asymptotic 
power law. The theory is easily extended to the aging 
regime where these secondary relaxations can dominate 
the rearranging motions. 

In RFOT theory, above the glass transition temper- 
ature, the entropic advantage of exploring phase space, 
manifested as a driving force for reconfiguration, is bal- 
anced by a mismatch energy at the interface between 
adjacent metastable states. For a flat interface in the 
deeply supercooled regime the mismatch energy can be 
described as a surface tension that can be estimated from 
the entropy cost of localizing a beadjH [TS], giving a sur- 
face tension a a = (3/4)/cBTr^" 2 ln[l/(c?| / 7re)] where <Il is 
the Lindemann length, the magnitude of particle fluc- 
tuations necessary to break up a solid structure, and 
is nearly universally a tenth of the inter-particle spac- 
ing, (<1l — O.lro). The free energy profile for reconfig- 
uration events resembles nucleation theory at first order 
transitions but is conceptually quite distinct. Following 
Stevenson-Schmalian-Wolynes (SSW)[]3] the free energy 
cost of an N particle cluster with surface area £ making 
a structural transition to a new metastable state may be 
written 

F(N,T,) = Ecro - Nk B Ts c - fc B Tlnfi(iV, £) - ^ Sf 

particles 

(i) 

A key element of the free energy profile is the shape en- 
tropy k B log Q(N, S) which accounts for the number of 
distinct ways to construct a cluster of N particles hav- 
ing surface area S. At one extreme are compact, nearly 
spherical objects with shape entropy close to zero. While 
objects such as percolation clusters or stringy chains have 
surface area and shape entropy both of which grow lin- 
early with N. The last term of equation [T] accounts for 
the inherent spatial fluctuations in the disordered glassy 
system that give fluctuations in the driving force. We 
presently ignore local fluctuations in the surface mis- 
match free energy, but their inclusion would not qual- 
itatively alter the results [13l [TBI El- We simplify by 
assuming uncorrelated disorder, so each particle joining 
the reconfiguration event is given a random energy, 5f, 
drawn from a distribution of width 8f. The r.m.s. mag- 
nitude of the driving force fluctuations above T g follows 
from the configurational heat capacity through the rela- 
tion 5f w Ty/AC p k B , a result expected for large enough 
regions. We will assume no correlations for simplicity, 
but they can be included. 



For nearly spherical reconfiguring regions forming com- 
pact clusters the shape entropy is very small and the mis- 
match free energy is (To47r(3A r /(47rpo)) e ^ 3 with 9 = 2 if 
fluctuations are small. In disordered systems the mis- 
match free energies grow with exponent 9 generally less 
than 2 reflecting preferred growth in regions of favor- 
able energetics and the large number of metastable states 
which can wet the interface and reduce the effective sur- 
face tension. A renormalization group treatment of the 
wetting effect [TS] suggests that 9 = 3/2 in the vicinity 
of an ideal glass transition. Incomplete wetting giving 
strictly 9 — 2 only asymptotically would not change the 
numerics of the present theory much. Whether com- 
plete wetting occurs for supercooled liquids under lab- 
oratory conditions is still debated[T71 [T5], [JJ5]. The free 
energy profile describing reconfiguration events restricted 
to compact clusters becomes, then, F compac t(N) — 
cr 47r(3iV/(47rpo)) e/3 - NTs c . The minimum number of 
particles participating in a reconfiguration event is de- 
termined by finding where the free energy profile crosses 
zero. For 9 — 3/2 the activation free energy barrier is in- 
versely proportional to the configurational entropy, lead- 
ing to the Adam-Gibbs [2"U] relation for the most probable 
relaxation time F$/k B T ~ lnr Q /ro ~ s^ 1 . Adding fluc- 
tuations to the profile of compact reconfiguration events 
yields an approximate Gaussian distribution of barriers 
with width scaling as V N^Sf. Xia and Wolvnes|12j. 
and more explicitly Bhattacharyya et al. |21j. have shown 
that, with the inclusion of facilitation effects, the result- 
ing barrier distribution accounts for the stretching of the 
main relaxation process and yields good estimates for 
how the stretching exponent varies with liquid fragility. 

Restricting the reconfiguration events to stringy clus- 
ters (using percolation clusters gives very similar results) , 
gives a free energy profile linear in the number of parti- 
cles reconfigured, save for the minimum cost Fj n to begin 
to reconfigure a region: 

F stlias (N) = -NT( Sc - S f in s) + F iD . (2) 

The critical "entropy" is given by Tsf nns — Vi n t(z — 
2) - k B T\n(z - 5) w 1.13fc s T. This is the difference be- 
tween the surface energy written in terms of the coordi- 
nation number of the random close packed lattice z ~ 12 
(v int is the surface tension per nearest neighbor) and the 
shape entropy including excluded volume effects [32]. If a 
bead can individually reconfigure then the cost to begin 
to reconfigure is Fi n = zvi nt — Ts c ~ 2.5 — 2.9k B T. If 
two must be moved then F m ss 6.1k B T. The continuous 
form of the surface mismatch energy gives a somewhat 
higher value when applied to these small regions, giving 

^continuous = ^4^3/(4^)6/3 _ Tf ,^ _ \Q^k B T for & 

one particle reconfiguration. The remarkably simple free 
energy profile of equation [2] monotonically increases, so 
that below T str ing (defined by s c (r str ing) = sf ring ) recon- 
figuration via pure stringy objects is impossible. Above 
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FIG. 1: The two dimensional free energy profile describing cooperative relaxation as a function of the number of compact 
particles JV C and the number of stringy particles Nf in the reconfiguring region. The transition state separating the unrecon- 
figured state (N c — Nf = 0) with the final stable state (colored yellow in the figure) determines the free energy barrier to 
reconfiguration. The upper left panel gives the profile (s c = l.Ofcs) in the absence of fluctuations while the others demonstrate 
three possible realizations of fluctuations for a relatively strong liquid, a liquid having ACp « Iks per bead. The fluctua- 
tions for the two dimensional profile are implemented as cumulative sums of local fluctuations in both iV c and Nf. For the 
situation described in the bottom right panel compact reconfiguration is required to overcome the free energy barrier. For the 
two other realizations the fluctuations are such that stable (yellow) regions exist along the vertical axis, meaning string-like 
reconfiguration is possible. These stringy clusters, stabilized by fluctuations, account for the secondary relaxation. 



^string the same process can occur with very small free 
energy barrier, having only to overcome F- ln . Thus T str i n g 
signals the crossover from dynamics dominated by acti- 
vated events to dynamics dominated by non-activated 
processes. Interestingly, this crossover is mathemati- 
cally analogous to the Hagedorn transition of particle 
theory [25] • The predicted constant value of s^ trin s i s con- 
firmed experimentally [14 . In contrast to the situation for 
compact reconfiguration, driving force fluctuations dra- 
matically alter the picture of stringy relaxation. With 
driving force fluctuations a lucky sequence of fluctuations 



can easily push the nominally linearly increasing free en- 
ergy profile to cross zero, so strings can be active below 
T c . 

SSWQ3] introduced a crude model to estimate the 
shape entropy and surface area of the range of shapes 
between the extremes. This "fuzzy sphere" model, con- 
sists of a compact center of N c particles with a stringy 
halo of Nf particles. With this interpolative model T g 
the preferred shape of a reconfiguring region is largely 
compact, while the relevant regions become more rami- 
fied close to T c . Fluctuations, aggregated cumulatively in 
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FIG. 2: Probability distribution of free energy barriers 
governing relaxation events in supercooled liquids. Differ- 
ent curves represent different temperatures, increasing from 
near the glass transition temperature to near the dynami- 
cal crossover temperature. The arrows indicate the typi- 
cal relaxation time predicted from the fuzzy sphere model 
without fluctuations. The dashed line gives the distribution 
of free energy barriers for a liquid just above the dynami- 
cal crossover temperature where primary relaxations disap- 
pear leaving only the secondary process. The top panel cor- 
responds to a rather strong liquid with small fluctuations, 
ACp w lfcp per bead — a material similar to Ge02- The 
bottom panel corresponds to a fragile liquid with larger fluc- 
tuations, ACp = 3fcp per bead — a material similar to ortho- 
terphenyl. Secondary relaxations, i.e. relaxation events us- 
ing string-like rearranging regions, increase in prominence as 
the temperature is increased, becoming the dominant pro- 
cess near the dynamical crossover temperature where the 
two peaks merge. The units of the activation energy are 
given in fcgT, which assumes a mismatch penalty primar- 
ily entropic in nature, ao ~ ksT . An energetic mismatch 
penalty, ao ~ feTx, would lead to Arrhenius behavior for 
the secondary relaxation process, as the distribution is peaked 
around the minimum free energy to initiate a stringy recon- 
figuration, _Fi n . In this calculation we have used the contin- 
uous approximation of F[ n shown in the text. The facilita- 
tion phenomenon, as described by Xia and Wolvnes|I2] and 
by Bhattacharya et al.[3T] but not accounted for here, would 
shift weight from the largest free energies to the center of the 
primary peak, raising the overall height of the primary peak 
relative to the secondary peak. 

the compact core yielding a variance proportional to N c , 
and cumulatively in the stringy halo yielding a variance 
proportional to Nt, modify the SSW two dimensional 
free energy landscape. The local free energy plots for 
several realizations of such accumulated fluctuations, as- 
sumed to have Gaussian statistics, are given in figure [I] 
For some realizations compact reconfiguration will still 
be required to overcome the free energy barrier, but for 
other realizations the fluctuations are such that the free 
energy cost for reconfiguration crosses zero along the Nt 
axis (shown in the figure in yellow) so the region is able 
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FIG. 3: Distribution of free energy barriers for a strong liquid 
(ACp ss Iks per bead) separated into the contribution from 
secondary relaxations (red curves), corresponding to string- 
like reconfiguration as in panels b and c of figure [I] and pri- 
mary relaxations (black curves), corresponding to compact 
reconfigurations. The full distributions are given for compar- 
ison (blue curves). The separation of the curves makes clear 
that as the dynamical crossover temperature is approached 
the primary relaxation becomes subordinate to the secondary 
relaxation. 



to relax via a string- like reconfiguration event. These 
stringy rearranging clusters, stabilized by disorder, we 
argue, are key contributors to the secondary relaxation 
process. The statistics of the stable reconfiguration paths 
with lowest free energy barrier are summarized in figure 
[2] with barrier distributions for two different values of 
Sf corresponding to a strong and a fragile liquid. We 
can dissagregate these distributions into the parts due, 
separately, to the compact events ("primary") and to 
the string-like fluctuation induced events ("secondary"). 
These distributions are shown in figures [3] and [4] At 
temperatures near T g compact primary relaxations dom- 
inate reconfiguration, but as the temperature increases, 
fluctuations are able to stabilize string-like reconfigura- 
tion more easily and the secondary relaxations increases 
in prominence. At the same time, with increasing tem- 
perature, the primary relaxation peak shifts to lower free 
energy barriers, making it difficult to distinguish between 
the two contributions. As T c is approached and crossed, 
the primary and secondary peaks merge with string-like 
reconfiguration clusters becoming the dominant mode of 
relaxation. For fragile liquids, i.e. liquids with larger con- 
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FIG. 4: The corresponding results to those of figure [3] but for 
a fragile liquid, one with ACp ~ 3kp per bead. 
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FIG. 5: Distribution of free energy barriers for the secondary 
relaxation process from the statistical sampling of the fuzzy 
sphere model with fluctuations. The data correspond to a 
strong liquid (ACp ~ 1) and show that at higher temper- 
atures (larger configurational entropies) the distribution de- 
cays more slowly. This leads to wider activation energy dis- 
tributions, matching the expectations of the analytical calcu- 
lations. 



figurational entropy fluctuations |10j . the secondary relax- 
ation peak is generally more important than for strong 
liquids. Both peaks are broader and begin to merge at 
lower temperatures for the more fragile liquids. Because 
facilitation effects are not explicitly accounted for (these 
would mostly affect the higher barriers) it is not easy to 
directly compare these predicted distributions with quan- 
titative precision to experiment. In addition, the number 
of reconfiguring particles in the two peaks is different, so 
their contributions to the measured amplitudes are differ- 
ent as well. Nevertheless, the predicted magnitude of the 
secondary relaxation peak, as compared to the primary 



peak, seems to be somewhat larger than experiments ap- 
parently show. This disparity is more pronounced for 
fragile materials. The assumptions in the fuzzy sphere 
model, and especially the assumption of uncorrclatcd dis- 
order, apparently overestimate the influence of the fluc- 
tuations, which are probably (anti-) correlated for the 
most fragile systems. 

The barrier distribution for reconfiguration events that 
take an ideal stringy form can be explicitly calculated. A 
similar analysis to the string case can be applied to recon- 
figuration via percolation clusters but with somewhat dif- 
ferent numerical constants in the relation of the free en- 
ergy profile to N. This analytic calculation resembles that 
of Plotkin and Wolynes for the "buffing" of protein fold- 
ing energy landscapes [24]. The key to the calculation of 
r^i* 1 *), the distribution of free energy barriers for events 
with less than N max displaced atoms, lies in mapping the 
problem onto a random walk, or a diffusion process in free 
energy space. Going to the limit of continuous number of 
particles we may write a stochastic differential equation 
for the free energy profile dF/dN = dF stl i n g/dN + Sf. 
The principal quantity to compute is G(N, F*; F; n ), the 
probability that a reconfiguration path of N particles has 
free energy F if the cost to initiate the reconfiguration 
event is F m . The evolution of G, that follows from the 
stochastic profile, is described by a diffusion equation 
with drift subject to absorbing boundary conditions at 
both F = and F = F*. These boundary conditions 
permit the calculation of the distribution of free energy 
barriers by keeping track of the maximum excursion of 
the random walk. 



dG 
dN 



dG 
dF 



V— 

2 J dF 2 ' 



(3) 



The slope of the mean free energy profile (j> — T(s^ trmg — 
s c ) depends simply on the proximity to the string transi- 
tion. is a string tension reflecting the free energy cost 
of lengthening a string. The probability density for the 
maximum excursion of F is then 



T(Ft) 



d 
dFi 



N„ 



dN 



2 



dG 
dF 



(4) 



F=0/ 0<F in <Ft 



The average (-)o<F ixl <Ft is present to integrate over the 
fluctuations in the free energy cost of initiating a string — 
capturing the statistics of the smallest possible activation 
barriers. The derivative with respect to F* converts from 
the cumulative probability that the free energy barrier 
is below F$ to the probability the free energy barrier is 
between F i and F* + dF* . G can be calculated explicitly 
by solving the diffusion equation using the method of 
images. The result may be represented in closed form in 
terms of the Jacobi theta function, however we leave the 
sum explicit to more easily examine the asymptotics 
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G = 



(F-F in -4>iV/2) 



y/2n6pN 



E 



2Sf 2 N 



(2nF$+F + F in ) 2 
2Sf 2 N 



(5) 



In the integral of equation [4] the cutoff N max reflects 
the maximum size to which a stringy reconfiguration 
event would typically grow before compact rearrange- 
ments dominate. We estimate this maximum length as 
JVmaj ps F$/(j), since certainly by that length the most 
important reconfiguration events would be compact. We 
can simplify by smoothing the cutoff so that the finite 
integral J^'^dN- is replaced by J^°dN exp(—N/N max )-. 
The smoothed distribution of free energy barriers follows 
directly from equations [4] and [5] 




exp 



sp 



-F in q 



exp(Fi n g) - 1 
exp(Ftq) - 



where, q = 




(6) 



0<F in <Ft 



The result involves nothing more complicated then expo- 
nentials and error functions. 

The total magnitude of the secondary relaxation peak 
is estimated by calculating the probability that fluctua- 
tions can stabilize a stringy reconfiguration for any size 
barrier less than F$. Integrating T over F$ yields 



W exp 



2<f>(F b 



(7) 



\P increases with temperature as the dynamical crossover 
at Turing is approached, a trend that is validated 
experimentally [25 . At the crossover temperature and 
above, this secondary relaxation becomes the only re- 
maining mode of activated relaxation. The sharp transi- 
tion from activated to non-activated motions at T c that 
is predicted by the non-fluctuating RFOT theory, as well 
as by mode coupling theory [57] and the mean field 
theory of supercooled liquids [25] [30], is therefore 
smoothed out by the string-like activated events made 
possible by fluctuations and exhibits no divergent crit- 
ical behavior at s^ tring [2T]. In this temperature regime 
the secondary beta relaxations of mode coupling the- 
ory would be present and overlap in frequency with the 
string-like activated secondary relaxations, perhaps mak- 
ing the differentiation of the processes difficult. 

r can be approximated by the Gumbel extreme value 
distribution function |31j. For F* > F[ n the barrier dis- 
tribution decays exponentially as 



r(F* > F in ) ~ exp(-F^) w exp 



_ 2 FH\ 

sp J 



(8) 



The results agree with the sampled distribution of barri- 
ers for the string-like reconfiguration events alone (shown 



. Using the fact that r = tq exp(F* /hsT) 
gives a power law distribution of relaxation 
~ t 1 where 7 pa 2(s s c 



in figure 
equation 

times P(t) ~ t~t where 7 ps 2(sf ring - s c )/AC P + 1. 
Well above T g the high barrier side of the secondary 
relaxation blends in with the primary relaxation peak. 
Thus the secondary relaxation from ramified reconfigu- 
ration events often appears as only a "wing" on the main 
distribution [521155] , 

In the aging glass, the picture of secondary relaxation 
is slightly modified. If the liquid fell out of equilibrium 
at Tf then the frozen-in structure has an average excess 
energy per particle e(T» = e(T K ) + dTAC P (T). At 
temperatures T < Tf a region of the liquid undergoing re- 
configuration would relax to a structure with average en- 
ergy e(T) < e(Tf). Thus the driving force for reconfigura- 
tion gains an energetic contribution and the configuration 
entropy in equation [I] is replaced by Ts c —> (Ts c + Ae) 
where Ae = e(T / )-e(T) = ffi dT'AC P (T'). Lubchenko 
and Wolvnes [TT] have shown that this additional driving 
force results in a change in slope of the typical relax- 
ation time as a function of temperature, and a transition 
to nearly Arrhenius behavior as the system falls out of 
equilibrium. Correspondingly, falling out of equilibrium 
causes the string tension to be reduced by an amount Ae, 
giving <j> — Turing _ rp s ^ _ ^ e anc j ma ]j m g the system 

appear closer to the dynamical crossover than an equi- 
librated system at the same temperature. Furthermore, 
the driving force fluctuations are frozen in as the aging 
glass falls out of equilibrium becoming largely indepen- 
dent of temperature. These changes broaden and flatten 
the barrier distribution as the temperature is lowered, 
giving 



T(F* >F-)~ exp ( 2 FHTs?^-T Sc -Ae) \ 

for large F*. The secondary relaxation strength in the 
aging regime becomes 



$ ps exp 



2 ( Ts st ring _ _ Ae ) 



TU B AC P (T f ) 



(F in - (Ts 



Ts c - Ae)) 



(10) 



In the limit T — > the distribution of barriers be- 
comes largely independent of temperature with T(F$ > 
F in ) ~ exp(-aFi) and a ps ({z - 2)v mt (T f ) - 
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Ae)/(T/fcsACp(T/)). For a broad enough distribution 
of barriers the dielectric absorption spectrum is deter- 
mined through the simple relation, e"(uj) ~ P(F$ — 



-fc^Tln w/wq) 



and becomes flat for low temper- 



atures, resembling the so called constant loss spectrum. 
In a rejuvenating glass, an aged system that is heated 
to a temperature above Tf , the energetic contribution to 
the driving force is negative, Ae < 0. In this situation 
the system appears as if it is be further from the dynam- 
ical crossover temperature than an equilibrated system 
at the same temperature and secondary relaxations are 
relatively suppressed. 

Nearly all glass forming liquids display, dynamical 
modes seemingly distinct from the primary alpha relax- 
ations. We have shown that by adding fluctuations to the 
existing structure of random first order transition theory 
a tail develops on the low free energy side of the acti- 
vation barrier distribution which shares many of the ob- 
served features of the secondary relaxations. The relax- 
ation process responsible for the tail differs from the pri- 
mary relaxation mechanism in the geometry of the region 
undergoing cooperative reconfiguration. While primary 
relaxation takes place through activated events involv- 
ing compact regions, secondary relaxation is governed by 
more ramified string-like, or percolation-like clusters of 
particles. While the existence of secondary relaxation 
is nearly universal, the relevant motions are of shorter 
length scales than those for primary relaxation, allow- 
ing additional material dependent effects and, perhaps, 
less universal quantitative description than for the main 
relaxation. The present theory, however, suggests a uni- 
versal mechanism for secondary relaxation. The theory 
points out some general trends about the way these relax- 
ations vary with temperature and substance which con- 
form to observation. 
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